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Abstract 

A kinetic consideration of an axisymmetric equilibrium with vanishing 
electric field near the magnetic axis shows that V/ should not vanish on 
axis within the framework of Vlasov theory while it can either vanish or 
not in the framework of both a drift kinetic and a gyrokinetic theories (/ 
is either the pertinent particle or the guiding center distribution function). 
This different behavior, relating to the reduction of phase space which leads 
to the loss of a Vlasov constant of motion, may result in the construction 
of different currents in the reduced phase space than the Vlasov ones. This 
conclusion is indicative of some limitation on the implications of reduced 
kinetic theories in particular as concerns the physics of energetic particles in 
the central region of magnetically confined plasmas. 
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There are important phenomena in plasma physics as microinstabilities 
which can not be described in the framework of a macroscopic theory as mag- 
netohydrodynamics but require the employment of kinetic theory. For the 
high temperature fusion plasmas an appropriate approximate kinetic theory 
is based on the Vlasov equation which the Boltzmann/Fokker-Planck equa- 
tion reduces to when the collision term is neglected. To solve self consistently 
the set of the Maxwell- Vlasov equations, however, is a tough problem related 
to the fact that the complete set of constants of motion is missing. For ex- 
ample, in the case of an axisymmetric equilibrium only the energy and the 
canonical momentum conjugate to the toroidal coordinate are known con- 
served quantities out of the four potential constants of motions. Because of 
the two missing constants of motion it is not possible to construct equilibria 
with macroscopic poloidal velocities, although phenomenologically sheared 
poloidal velocities play an important role in the transition to improved con- 
finement regimes of magnetically confined plasmas, e.g. the L-H transition. 
In addition, the temporary and probably future computational efficiency puts 
a limitation on numerical solutions of the Maxwell- Vlasov equations. For this 
reason approximate kinetic theories in a reduced phase space, as the drift ki- 
netic [lj-[6] and gyrokinetic ones [7]- [15] have been developed and applied to 
numerous simulations, e.g. on turbulence in connection with the creation of 
zonal flows. In both theories the reduced phase space is five dimensional with 
three spatial components associated with the guiding center position, R (or 
the gyrocenter position in the framework of gyrokinetic theory), and a veloc- 
ity component parallel to the magnetic field, v\\; also, the two compoments of 
the perpendicular particle velocity are approximated after a gyroangle aver- 
aging with the magnetic moment which is treated as an adiabatic invariant. 
A related underlying assumption for both reduced theories is that the ratio e 
of the gyroperiod to the macroscopic time scale is small. In the drift kinetic 
theory e is the same as the ratio of the gyroradius to macroscopic scale length 
while in the gyrokinetic theory small spatial variations are permitted but, for 
example, the amplitudes of the fluctuations to the background fields is equal 
to e [Hi]. It may be noted here that the reduced-phase-space kinetic theories 
are developed via expansions in e the convergence of which is not guaranteed. 

Because of the reduction of the phase space some information of the 
particle motion is missing. This gives rise to the question: is the missing 
information important? In the present note we address this question by 
making a comparison of Vlasov with drift kinetic and gyrokinetic theories 
near the magnetic axis of an axisymmetric magnetically confined plasma 
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with vanishing electric fields. Motivation was a previous study [TB] in which 
by considering this equilibrium in the framework of Vlasov theory we found 
the following new constant of motion: C\ = v z + / In where v<j, is the 
toroidal particle velocity, v z the velocity component parallel to the axis of 
symmetry and = I/rh^ the toroidal magnetic field near axis (r, <p, z are 
cylindrical coordinates). For the sake of notation simplicity and without loss 
of generality here we will consider only ions and employ convenient units by 
setting m = q = c = 1 where m and q are the ion mass and charge and c is 
the velocity of light. Since phenomenologically the density gradients vanish 
on the magnetic axis it was a surprising conclusion of Ref. [16] that V/ ^ 
must hold on axis, where / is the particle distribution function. The reason 
is that if one assumes V/ = thereon it is not possible to obtain one of the 
known constants of motion, i.e. the canonical toroidal momentum C2 = rv^. 
Note also that because of the absolute value of in C\, distribution functions 
depending on C\ and the energy C3 = l/2[v^+vl+v^] can not create currents 
parallel to the magnetic field. 

The form of C\ relates crucially to the toroidicity because in the straight 
case for which changes to z, z to y and r to x (where x, y, z are Carte- 
sian coordinates) for a ^-independent equilibrium with straight magnetic axis 
parallel to z and arbitrary cross sectional shape the respective constant of 
motion becomes C\ = v z . Thus, in this case distribution functions depending 
on C\ and the energy can produce parallel currents. It may also be noted 
that in the straight case for V/ 7^ on axis one can obtain the complete set 
of four constants of motion near axis: C\ = v z , C2 = v x — B z y, C3 = v y + B z x 
and C4 = l/2(vl + v y +v^) where B z is the "toroidal" magnetic field on axis 
[T6j . Therefore, distribution functions of the form /(C2, C3, C4) can lead to 
purely "poloidal" velocities near axis irrespective of the cross sectional shape 
in consistence with magnetohydrodynamics [TB] . Unlikely, toroidal magneto- 
hydrodynamic equilibria with purely poloidal velocities are not possible [T9] . 
Whether this conclusion of nonexistence persists in the framework of Vlasov 
theory is an open question relating to the two unknown constants of mo- 
tion. The above comparison shows that the Vlasov theory well distinguishes 
equilibria with circular and straight magnetic axes. 

In the present note we examine the same equilibrium near axis in the 
framework of drift kinetic and gyrokinetic theories on an individual basis. 
Though the drift kinetic equations of Ref. [3] and the gyrokinetic equations 
of Ref. will be employed we claim that the conclusions do not rely on the 
particular forms of the reduced kinetic equations. 
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Drift kinetic theory 



The drift kinetic theory established in Ref. j3] is based on the Little- 
john's Lagrangian for the guiding center motion [T7] extended to include the 
polarization drift in such a way that local conservation of energy is guaran- 
teed. The drift kinetic equation for the guiding center distribution function 
/(R, v\\, /i, t) (with ft = 0) acquires the form 

df . df 

The guiding center velocity, v, and the "acceleration" parallel to the magnetic 
filed, v», can be expressed in a concise form by introducing the modified 
potentials, A* and $*, and the respective modified electric and magnetic 
fields, E* and B*, as [I]: 



A* = A + ^||b + v B , (2) 

$* = $ + ^ + 1(^ + 4), ( 3 ) 

EXB 

= -gs-, (4) 

E * = -er-w (5) 

B* = VxA*, (6) 

where $ and A are the usual electromagnetic scalar and vector potentials 
and b = Ti/B. The quantities v and v\\ are then given by 

B* E*xb 

V = W 9 = V \\^ + — ( 7 ) 



E* B* 1 



-v 



E*, (8) 



where 



Sf = B*-b = B + U||b- Vxb + b VXv E . (9) 

Explicit expressions for v and are given by Eqs. (3.24) and (4.17) of Ref. 
[3]- Also, it is noted here that Eqs. (J7J) and (jSJ) have similar structure as the 
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respective gyrokinetic equations of Refs. P2] and [IT] [Eqs. (5.39) and (5.41) 
therein]. Since B* appears in the denominators of (JTj) and (jHJ) a singularity 
occurs for Bi — 0. For E = this singularity can be expressed by the 
critical parallel velocity v c = —Q/(h ■ Vxb), where Q is the gyrofrequency. 
Therefore, the theory is singular for large \v\\ \ at which v and v\\ diverge and 
consequently non-casual guiding center orbits occur and the guiding center 
conservation in phase space is violated [2J. It is the -dependence of A* [Eq. 
02])] that produces the singularity. In order to regularize the singularity v\\ in 
(J2j) can be replaced by an antisymmetric function g(z) with z = v\\/vq, where 
Vo is some constant velocity [2]-[l]. The nonregularized theory presented here 
for simplicity is obtained for g(z) = z. In the regularized theory g(z) ~ z 
should still hold for small \z\. For large \z\, however, g must stay finite such 
that with v 3> ^thermal one nas v o9(°°) < v c- A possible choice for g(z) is 
g(z) = tanh z. 

As in the Vlasov case [16] we consider the drift kinetic equation (JTJ) for 
an axisymmetric equilibrium with E = in the vicinity of the magnetic axis. 
Since on axis the magnetic field becomes purely toroidal and dependent only 
on r one readily calculates 



Vxb = Vxe„ 



2l 

r 



Bt = B, B* = Bej, + v\\ — , Vfi(r) = ^e r 

ii r ^ r 



and consequently 



(v\ u dB\ . . 

v = v fl = , l|e , + ---- e z , (10) 



U||=0. (11) 

As expected on axis the guiding center velocity consists of a component 
parallel to B and the curvature and grad-S drifts perpendicular to B and 
parallel to the axis of symmetry. Also, the "acceleration" vanishes because 
there is no parallel force and the drift kinetic equation ([]]) becomes 

vf^O. (12, 

Therefore, unlike the Vlasov description near axis the distribution function 
because of ( TTTT) can be any function of v\\ , which is a constant of motion, and 
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on axis can hold either V/ ^ or V/ = 0. Thus, irrespective of the value 
of V/ on axis one can either construct parallel currents or not by choosing / 
either symmetric or antisymmetric in v\\ . This is a significant difference from 
the Vlasov situation in which, if V/ = on axis the obtainable particle dis- 
tribution functions of the form f(Ci,Cs) can not produce parallel currents. 
This discrepancy clearly relates to missing C\ in the reduced phase space 
which results in an nontrivial loss of information for the particle motion. 
Note that near axis the overwhelming majority of the particles are passing 
and the parallel currents constructed in the framework of the drift kinetic 
theory may differ from the "actual" ones. In the context of the drift kinetic 
theory the two constants of motion, i.e. the energy fiB + l/2un and the 
canonical momentum r(A^ + vnb^) can be found from ([T]) everywhere by the 
method of characteristics. Thus, the complete set of constants of motions 
is obtained in the five dimensional phase space. These are recoverable on 
axis where v», r, A<j,(r) (and fi) are conserved individually even if V/ = 0. 
Unlikely, in the respective Vlasov case the toroidal angular momentum con- 
stant of motion, C*2, is missed when V/ = is assumed. Also, for straight 
^-independent equilibria one has on axis Vxb = Vxe 2 = 0, B* = Be z , 
v 9 = v\\e z , v\\ = and (JT]) is identically satisfied, implying that / can be any 
function of x, y, vn and fi. Therefore, unlike the Vlasov theory the depen- 
dence of / on v || is independent of toroidicity and therefore, regarding the 
formation of parallel currents, the drift kinetic theory can not distinguish 
equilibria of circular or straight magnetic axes. 

Gyrokinetic theory 

We will use the gyrokinetic equations of Ref. [7] which have been em- 
ployed to a variety of applications (see for example the recent Refs. [121 
H~3| IT4"]). Eq. flTJ remains identical in form where /(R, v\\, t) is now the 
gyrocenter distribution function for ions. The gyrocenter velocity and "ac- 
celeration" are given by 



v = v g = v\\b + 



B 



(V B + V VB + V c ) , 



(13) 




(14) 
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Here, B is the equilibrium magnetic field, b = B /S , 

5*1, = (Bo + t;||Vxb ) -b , 

$ stands for the perturbed gyroaveraged electrostatic potential, and the 
ExB-drift velocity Vg, the grad-£> drift velocity v V s, and the curvature 
drift velocity v c are given by 



Note that as in the case of drift kinetic theory a similar singularity occurs 
at 5q|| = 0. In numerical applications this singularity was "avoided" by 
approximating B£,< = Bq (see for example Refs. (T2J [13]). Consideration of 
the above equations for an axisymmetric equilibrium with E = near axis 
yields relations similar to (JTQJ) , (II ip and ( fT2l) . Therefore, the above found 
discrepancies of the drift kinetic theory with the Vlasov one persist in the 
framework of the gyrokinetic theory. The structure of the reduced kinetic 
equations in conjunction with the symmetry of the equilibrium considered 
clearly indicate that this conclusion is independent of the particular drift 
kinetic or gyrokinetic equations. 

In conclusion, first a singularity which occurs in both drift kinetic and 
gyrokinetic theories for large parallel particle velocities is usually eliminated 
in the literature by a rough approximation. Second, a comparison of the 
Vlasov equation with either the drift kinetic or the gyrokinetic equation near 
the magnetic axis of an axisymmetric equilibrium with vanishing electric field 
implies different properties of V/ and, unlike Vlasov, non distinguishing of 
equilibria with straight and circular magnetic axes in connection with the 
formation of parallel currents. This relates to the loss of a Vlasov constant 
of motion in the reduced phase space. Consequently, different drift kinetic 
or gyrokinetic parallel currents may be created than the Vlasov ones. This 
indicates that the reduction of the phase space, even if made rigorously so 
that local conservation laws and Liouvillean invariance of the volume element 
is guarantied, is associated with the loss of nontrivial physics which could 
put certain limits on the validity of the conclusions from drift kinetic or 
gyrokinetic simulations. 



v$xva 



(15) 





(16) 



(17) 



7 



Acknowledgements 



Part of this work was conducted during a visit of G.N.T. to the Max-Planck- 
Institut fur Plasmaphysik, Garching. The hospitality of that Institute is 
greatly appreciated. 

This work was performed within the participation of the University of 
Ioannina in the Association Euratom-Hellenic Republic, which is supported 
in part by the European Union and by the General Secretariat of Research 
and Technology of Greece. The views and opinions expressed herein do not 
necessarily reflect those of the European Commission. 



8 



References 

D. Pfirsch, Z. Naturforsh. 39a, 1 (1984). 

D. Correa Restrepo and H. K. Wimmel, Physica Scripta 32, 552 (1985). 

D. Correa Restrepo, D. Pfirsch and H. K. Wimmel, Physica 136A, 453 
(1986). 

D. Pfirsch, P. J. Morrison, Phys. Fluids B 3, 271 (1991). 

G. N. Throumoulopoulos and D. Pfirsch, Phys. Rev. E 49, 3290 (1994). 

H. Weitzner, Phys. Plasmas 2, 3595 1995. 
T. S. Hahm, Phys. Fluids 31, 2670 (1988). 
H. Sugama, Phys. Plasmas, 7 466 (2000). 
A. Brizard, Phys. Plasmas, 7 4816 (2000). 

D. Pfirsch and D. Correa-Restrepo, J. Plasma Phys., 70 719 (2004). 

D. Correa-Restrepo and D. Pfirsch, J. Plasma Phys., 70 757 (2004). 

X. Lapillonne, B. F. McMillan, T. Gorier, S. Brunner, T. Dannert, F. 
Jenko, F. Merz, and L. Villard, Phys. Plasmas 17, 112321 (2010). 

Feng Liu, Z. Lin, J. Q. Dong, and K. J. Zhao, Phys. Plasmas 17, 112318 
(2010). 

H. S. Zhang, Z. Lin, H. Holod, X. Wang, Y. Xiao, and W. L. Zhang, 
Phys. Plasmas 17, 112505 (2010). 

Lu Wang and T. S. Hahm, Phys. Plasmas 17, 082304 (2010). 

H. Tasso and G. N. Throumoulopoulos, J. Phys. A: Mathem. Theor. 40, 
F631 (2007). 

R. G. Littlejohn, J. Plasma Phys. 29, 111 (1983). 

G. N. Throumoulopoulos and G. Pantis, Plasma Phys. Controlled Fusion 
38, 1817 (1996). 

[19] G. N. Throumoulopoulos, H. Weitzner and H. Tasso, Phys. Plasmas 13, 
122501 (2006). 



9 



